

  fna = "sgp_2010_cam5_v2.nc" 
  fnb = "sgp_20100401_shaocheng.nc" 

  fla = addfile(fna,"w") 
  flb = addfile(fnb,"r") 

  bdate = flb->bdate
  bdate = (/ 20100301 /) 
  fla->bdate=bdate 

  lat = 36.605 

  lat!0 = "lat" 
  lat&lat = 36.605

  lon = -97.485

  lon!0 = "lon" 
  lon&lon = -97.485
 
  fla->lon=lon
  fla->lat=lat
  time = fla->time 

  lev = flb->lev
  fla->lev=lev

  nt = dimsizes(time) 
 
  tsec = new(nt,"double") 
  tsec!0 = "time" 
  tsec&time = time
  tsec@long_name = "time" ;
  tsec@units = "seconds" ;

  tsec = (/ fspan(0.,3600.*2207,2208)   /) 

  fla->tsec=tsec 


 
